Estimating the conditions for polariton condensation in organic thin-film microcavities 
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We examine the possibility of observing Bose condensation of a confined two-dimensional polari- 
ton gas in an organic quantum well. We deduce a suitable parameterization of a model Hamiltonian 
based upon the cavity geometry, the biexciton binding energy, and similar spectroscopic and struc- 
tural data. By converting the sum-over-states to a semiclassical integration over d-dimensional phase 
space, we show that while an ideal 2-D Bose gas will not undergo condensation, an interacting gas 
with the Bogoliubov dispersion H{p) ~ sp close to p = will undergo Bose condensation at a given 
critical density and temperature. We show that Tc/^/fh- is sensitive to both the cavity geometry and 
to the biexciton binding energy. In particular, for strongly bound biexcitons, the non-linear inter- 
action term appearing in the Gross-Pitaevskii equation becomes negative and the resulting ground 
state will be a localized soliton state rather than a delocalized Bose condensate. 



Over the past few years, considerable progress has been 
achieved in creating Bose Einstein condensates (BEC) 
of weakly interacting particles. The most recent break- 
through is where BEC of photons was realized by pump- 
ing an optical cavity filled with a dye and bounded by 
two reflective mirrors [l]. Within the optical cavity, the 
photons acquire an effective mass as determined by the 
cut-off frequency of the cavity that can be 6-7 orders 
of magnitude less than mass of an electron. Depending 
upon the density, this allows for a BEC transition tem- 
perature that can approach room temperature. Polari- 
tons are also ultra-light quasiparticles that are known to 
condense in systems composed of a semiconducting quan- 
tum well sandwiched between two reflective mirrors [5*- 
[S]. In this case, however, the polaritons act as hard-core 
Bosons and scattering at high density allows for a rapid 
thermalization of the gas. Recent experiments by Kam- 
mann et al. report on the cross-over between photon and 
polariton BEC within the same microcavity device. [5] . 

In this paper, we consider the formation of a polariton 
BEC in an organic quantum well consisting of a layer of a 
semiconducting polyacene molecular crystal sandwiched 
between two reflective mirrors. In contrast to inorganic 
quantum wells, organic quantum wells can be fabricated 
on the bench-top without the need for ultraclean environ- 
ments. Recent experiments by Forrest's group indicate 
that exciton polaritons in organic quantum wells can ex- 
hibit a lasing transition at room temperature [7*, and 
our recent work indicated that parameteric amplification 
and the transition to superfluid states should be possi- 
ble in these systems even in the presence of some disor- 
der within the crystal lattice itself |9j. Furthermore, as 
we shall discuss, the fact that molecular crystals may be 
glassy or have finite-sized microcrystaline domains breaks 
translation symmetry and may enhance the local density 
of polaritons thereby resulting in a higher BEC transi- 
tion temperature. We consider the nature of the po- 



lariton/polariton scattering interaction by relating the 
biexciton binding energy to the S'-wave scattering length 
for excitons in the cavity. Finally, taking molecular-level 
data as input to our theory, we estimate the transition 
temperature and critical density needed to achieve polari- 
ton BEC in a quasi-two dimensional molecular crystal. 



Polaritons are quasiparticles formed by strong coupling 
between a photon field and excitations within a material. 
Within the Frenkcl exciton model, wc can write a model 
Hamiltonian for this assuming that excitations are local 
to molecular sites, that excitons can hop to other sites 
via dipole-dipole coupling of their transition moments, 
and that the exciton/photon interaction in mediated 
via single quanta exchanges. Since polaritons are very 
lightweight quasi-particles, it is best to work in the long- 
wave, continuum limit rather than within a molecular- 
site representation more suitable for Frenkel excitons. 
However, to do so, we need to determine the effective 
mass of both the excitons and cavity photons. The co- 
herent motion of Frenkel excitons between sites is deter- 
mined by coupling between transition dipoles on different 
molecular sites, J{r). Since the dipole/dipole interaction 
is short ranged and centrosymmetric J{r) — J{—r), the 
coherent motion of the exciton in the Wigner representa- 
tion reduces to 1/F(r,p) = {p/rriex) ■'^rW{r,p) where the 
effective mass of the exciton is given by summing over all 
sites = J^n J(.^n)rn/^, where A is the area of the 
2D slab[10j. This allows us to readily compute the exci- 
ton effective mass given an instantaneous arrangement of 
the molecules within the 2D quantum well. Using the ex- 
perimental crystallographic arrangement of anthracene, a 
So — >■ 5*1 vertical transition moment of fii = 1.1901e bohr 
and a dielectric constant of e = 2.99 [11 we obtain an 
effective mass of — 6.83me which is close to the 
value of lOme typically reported for polyacene molecular 
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FIG. 1: Anthracene lattice layer from MD simulation. The 
repulsive core radius, a, and biexciton binding region are in- 
dicated by the concentric circles. The sketch to the right in- 
dicates the effective square-well potential for exciton/exciton 
scattering. 



crystals. |12j Thus, we write the exciton dispersion as 



2m„ 



where hio'^^ is the exciton transition energy. For an- 
thracene, we take this to be 3.1117 eV corresponding to 
the experimental vertical absorption energy [T3j. Within 
the cayity, the photon dispersion is given by 



where L is spacing between the cavity mirrors and r/ = ^/e 
is the refractive index. If we expand this for small values 
of k, one obtains the parabolic dispersion 



2m 



(1) 



where A = 2Trhc/Lri is the cut-off energy for the cavity 
and rrieff — 27:1] /cL is the effective photon mass. Tak- 
ing L = 250nm, -q = 2.99[IT] one obtains A = 2.87eV 
and rrieff/me = 1.67 x 10~^. These conditions put the 
cavity cutoff in resonance with the exciton transition en- 
ergy, thereby producing the strongest coupling between 
the photon field and the excitons. 

In the limit of low density, the combined exciton/cavity 
system is described by the energy functional 

f dr¥{r).H,.^{r)+g, f dr{^l;l{r))\Mr))\ {2) 



where 'l'^(r) = (-0^2. (r), ■^^^^^(r)) is a two-component field 
containing both exciton and photon contributions and 



H„ = h 



n/2 OJcav 



V{r), 



(3) 



describes the mixing between the exciton and photon 
fields. Note, that we have included a trapping potential, 
V{r) that may include any external effect that breaks 
translational symmetry. 

Polariton/Polariton kinematic interactions were stud- 
ied in detail by Zoubi and La Rocca [H] using the 
Agranovich-Toshich transformation |15j which converts 
the Frenkel excitons to Bosons, thereby allowing the exci- 
ton/exciton scattering interaction to be derived. In short, 
they show that in the long-wave limit and in the case 
where the LP dispersion is nearly parabolic, the scat- 
tering T-matrix between polaritons is identical to that 
scattering from a square well with radius a. With this in 
mind, we estimate polariton/polariton coupling by as- 
suming that this is due to exciton/exciton scattering. 
However, here we shall construct an effective square- well 
potential for exciton/exciton scattering based upon the 
molecular geometry and available spectroscopic data. 

The non-linear coupling coefficient, — A'Kha/mex 
is determined by the low-energy S-wave exciton/exciton 
scattering length, a. To provide a robust estimation of 
this term, based upon molecular and spectroscopic prop- 
erties, we assume that stable bi-exciton states exist as a 
single bound state in an effective square well potential 
with binding energy Eb and that the range is determined 
by the nearest neighbor spacing in the molecular lattice, 
h as sketched in Fig. [T] Since excitons are Fermions in 
that we must restrict their population to at most a sin- 
gle exciton per local anthracene site, we need to include a 
hard-core contact radius, a, which we take to be the van 
der Waals radius of anthracene in the molecular crystal 
plane. From quantum scattering theory we arrive at a 
simple expression for exciton/exciton S-wave scattering 
length 



tan((6 — a)A;,] 



a + (5 - a) 1 



[b - a)Xb 



where A;, = \/2mex£b/fi- Since a may be either positive 
or negative depending upon the the bi-exciton binding 
energy, it is possible that different organic media within 
the cavity would give rise to either bright (a < 0) or 
dark soliton {a > 0) ground states. Taking b = 2a, the 
crossover between bright and dark soliton states would 
occur when the de Broglie wavelength of bound biexciton 
state: A;, = h/yy2mex£b ~ a/\/3. For systems such as 
anthracene with Sb ~ lOOcm^^ and taking a to be 3.5A, 
Afc/a w 2.2 indicating that we would expect that the 
hard-core interaction to dominate, producing a > for 
most polyacene systems. It is necessary to point out, 
however, that the presence of an attractive long-range 
exciton/exciton interaction does not immediately imply 
that there will be a bound bi-exciton state. It is readily 
shown that in 1 and 2 dimensions, an arbitrarily shallow 
square well will always support at least 1 bound state; 
however, in 3 dimensions, this is not the case. 



3 



Polaritons are formed by transforming Ho into a diag- 
onal representation, yielding upper and lower branches 
with dispersion 



± 



where 51 is the Rabi frequency and 9k = 51/ (w™'" — ) is 
the mixing angle. If we ignore mixing between the upper 
and lower branches (due to exciton/exciton scattering), 
the ground state of the system is a stationary solution of 
the Gross-Pitaevskii (GP) equation [T7] . 



[hwLP + V{r) + g|0c(r)|') 0c(r) = n4>c{r), 



(4) 



where /i is the chemical potential. In this expression, 
the ipc is normalized to unity since we have included the 
number of excitations in system into the non-linearity, 
g = gx^o-^cav, where is the Hopfield coefficient for 
the exciton portion of the LP wave function evaluated at 
fc = and Ncav is the number of photons in the cavity. 

For an ensemble of Bosons, the average number of par- 
ticles in the gas is given by 



TV 



A 



1 - A 
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(5) 



where A — exp(/i/3) and we have pulled out the term rep- 
resenting the population of the ground state. The sum 
in Nex is over all excitation of the system. An ideal gas 
of bosons will not form a condensate unless the dimen- 
sionality of the system, d > 2. The reason for this is that 
density of states scales as u}{e)de oc e'^^'^~^de and con- 
sequently the sum over the excited states diverges and 
BEC is forbidden by translational symmetry in 1 and 2 
dimensions for a non-interacting system. 

In the experiments by Snoke and coworkers [U [18], 
translational symmetry is broken by putting a local strain 
on the cavity. It has also been argued that the finite size 
of the pumping laser is sufficient to break the 2D trans- 
lational symmetry. Likewise, in the photon condensation 
reported by Klaers et al., condensation occurred within 
a finite 2D cavity with parabolic mirrors [1]. For the case 
of a molecular crystal, translational symmetry can be 
broken by the presence of microcystaline domains within 
the sample itself. Within a given domain, the electronic 
transition moments are aligned along the crystallographic 
axes. The grain boundaries impose a hard-wall trapping 
potential. This should allow for a much higher density of 
polaritons since excitons created within a given domain 
will be confined to that domain rather than being free 
to diffuse through out the sample. However, as we show 
next, grain boundaries alone are insufficient to break 2D 
symmetry and one needs to consider the effect of the in- 
teractions. 

In order to estimate the transition temperature for a 
polariton gas in an organic microcavity, let us assume 



that the number of polaritons is high enough that we 
can ignore the kinetic energy term in the GP equa- 
tion and that V{r) is a circular well potential of ra- 
dius Ro corresponding to the typical size of an crys- 
talline microdomain. Within the well, we can write 
(/)^ = (l/(27ri?2))i/2 and thus, ^ = g/{2TTRl) (recah, that 
g already includes the number of quasi-particles in the 
system). Excitations of the condensate are then deter- 
mined by writing ip — (j)c + S^p and reintroduce this into 
the GP equation. Assuming the energy dispersion for 
the non-interacting system is quadratic in k the energy 
dispersion for the excitation is given by the Bogoliubov 
dispersion 



1/2 



For g sufficiently large, we can approximate the Bogoli- 
ubov spectrum as — s\p\ where s is the speed of sound 
s = {g/niLy^'^ and p — hk is the momentum. In this 
regime where the dispersion is linear in p, the conden- 
sate is behaving as a superfluid according to the Landau 
criterion. 

To evaluate Ng^ , let us write the classical Hamiltonian 
for a particle in a spherical potential well with radius Rg 
using H(p, r) = sp'^ and invoke the semi-classical approx- 
imation. For a = 1 we have the linear dispersion for the 
Bogoliubov excitation and for a = 2 we have the ordinary 
quadratic dispersion. Thus, we can convert the sum over 
states to a integral over phase space. 



p'^-^dp 



-^dr, 



(6) 



with Cd = 7r''/^/r(d/2 -|- 1). The resulting integral con- 
verges for all dimensions and positive values of U and d. 
In particular, for a 2D cavity and a = 1 with radius Ro, 



N, 



ex 



2c2Li2(A) 



(7) 



where A is the surface area of the microdomain and 
Li„ (x) is the polylog function. Eq. [7] allows us to predict 
the transition temperature for the formation of a polari- 
ton BEC using input strictly from molecular or spectro- 
scopic considerations. It follows that, 



N 



= 1 - 



for T <To and 



(8) 



(9) 



It is important to note that the cavity properties and 
number of cavity photons enters into this relation in two 
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ways. First in the density, p w X'^pcav and in the sound 
velocity, s. Fig. [2] shows the variation in the transition 
temperature (scaled by ^Jpcav) for various values of the 
biexciton binding energy and the cavity size. The vari- 
ation with cavity size reflects the exciton fraction of the 
polariton wavefunction. Close to resonance, the Hop- 
field coefficient, X^, is at its maximum value. It is in- 
teresting to note that as the biexciton binding increases, 
the transition temperature generally decreases. This is 
because the attractive part of the exciton/exciton inter- 
action decreases the effective exciton scattering length. 
However, for strongly bound bound bi-excitons, the scat- 
tering length can become negative resulting in an transi- 
tion from a BEC state to a bright-soliton state. In this 
case. To becomes imaginary and the transition tempera- 
ture can not be defined. For the model cavity at hand, 
this occurs when sb — 26.6meV. This is a large biexciton 
binding energy for organic molecular crystal systems. 

In conclusion, we predict that polariton BEC should 
be readily observable in an organic microcavity system 
using polyacene thin films. The arguments we present 
are based entirely up either molecular or spectroscopic 
parameters. Our analysis is based upon the notion that 
the non-linear interaction in the GP equation can be de- 
duced from the biexciton scattering which we treat in 
terms of a repulsive inner core surrounded an attractive 
square-well interaction. For strongly bound biexcitons, 
the exciton/exciton scattering length can become nega- 
tive and the resulting polariton ground state will be a 
bright-soliton rather than a Bose Einstein condensate. 
The results presented here are not limited to organic 
molecular crystals and apply equally to J-aggregate sys- 
tems, such as those recently studied by the Bulovic group 
at MIT [IS]. 
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FIG. 2: Scaled BEC transition temperature versus biexciton 
binding energy for various cavity sizes 
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